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INTRODUCTION 


This final report includes work performed on the "Study of the 
Mechanisms of Flux Pinning m Type-II Superconductors,” under NASA 
contract NAS8-24816. The work reported was performed during the first 
year of a project expected to be of two to three years duration During 
the course of this work, it became necessary to limit the program to 
the funds allocated for one year of effort. 

The term "flux pinning" refers to the interaction of quantized 
vortices of magnetic flux (fluxoids) m type- I I superconductors with 
defects in the superconducting sample This interaction restricts the 
motion of fluxoids in the sample and can pm them to the sites of de- 
fects. Since fluxoids are restricted in their motion, a gradient m 
the flux density results and a current density is established such that 
the pinning force equals the Lorentz force on the fluxoids 

In materials without defects, a current density cannot be estab- 
lished. Thus, defects are introduced intentionally to obtain materials 
with currents useful for practical purposes. Generally, defects are 
introduced m superconductive alloy wires m the process of reducing 
them from ingot to wire size. Superconductive mtermetallic compounds 
are produced generally at high reaction temperatures and normally con- 
tain many defects. Some degree of control over the defects introduced 
is possible in the metallurgical processing of these materials, but un- 
fortunately at present there is no method of determining which defects 
are important in establishing a specific magnitude or characteristic of 
the resulting current density. Present materials were developed 
empirically. 
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The objective of this work was to develop a method of determining 
the effect of various types of defects on the current density and thereby 
specifying ways of making superconductors with properties appropriate to 
specific applications. 

The fundamental property of defects used m this program was the 
change in the Gmzburg-Landau parameter H associated with each type of 
defect. Knowing this change m k and the distance over which the change 
occurs, the change in energy with distance, 1 e., the pinning force, can 
be determined 

The program consisted of three parts (l) development of a theo- 
retical model capable of describing the interaction of fluxoids with 
specific types of defects, ( 2 ) determination of the distribution of n 
m the superconductor using measured values of the specific heat 
capacity,* and ( 3 ) measurement of the distance between fluxoids 
(from which fluxoid forces can be determined) and determination of the 
underlying defect structure using electron microscopy. From (3) above, 
the relative numbers of defects of various types would be evaluated and 
from (2), the Ah associated with each defect would be determined 

Parts (l) and (2) above were completed successfully, but experimen- 
tal difficulties were encountered m part (3) that prevented us from 
obtaining the data required to treat the entire flux pinning problem 
coherently and verifying the theoretical predictions. 


The specific heat measurements were made at Stanford University by 
Mr Robert Zubeck and are being continued under the direction of 
Dr T H Geballe, These measurements, together with their interpre- 
tation, are expected to form the basis of Mr Zubeck's Ph D. dissertation. 
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SUMMARY 


Theoretical and experimental studies were made of the flux pinning 
mechanisms m type-II superconductors* A theoretical model based on 
the London theory was used to obtain trial functions for a variational 
calculation using the Gmzburg-Landau theory The microscopic proper- 
ties of fluxoids were calculated* The model was applied to the problem 
of flux pinning by (l) large defects having a h gradient, ( 2 ) small 
normal regions or voids, and (3) small regions having k values different 
from the bulk of the sample The forces arising from these pinning 
mechanisms were related to the critical transport current and were 
found to provide reasonable approximations to experimental data Item 
(3) above seems to be responsible for the peak effect m type-II super- 
conductors • 

Specific heat measurements were made on annealed and deformed pure 
niobium samples Changes m the normal- superconducting transition in- 
duced by straining the samples were analyzed, yielding the distribution 
of k and T Q throughout the sample The method developed should be val- 
uable for the general analysis of defects in superconducting materials 

Attempts were made to decorate fluxoids using the Bitter pattern 
technique and correlate the observed distributions with the underlying 
defect structures* Tight adhesion of the ferromagnetic powder to the 
sample prevented stripping the powder for observation m an electron 
microscope • 
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THEORETICAL MODEL 


Introduction 

Many properties of defect- free type- II superconductors have been 

1 2 

described by Abrikosov, Gorkov and others based on the -Ginzburg- 
3 

Landau theory of superconductivity. One of the major distinctions be- 

4 

tween the Ginzburg-Landau theory and the earlier London theory of 
superconductivity is the spatial dependence of the order parameter, ijr, 
m the Ginzburg-Landau theory and its constancy m the London theory 
This difference affects the predicted behavior of the density of super- 
conducting electrons m the vicinity of a boundary between normal and 
superconducting regions of the material. The London theory predicts a 
uniform condensation of electrons into the superconducting state while 
the Ginzburg-Landau theory predicts a condensation gradient of super- 
conducting electrons over a distance characterized by the coherence 
distance ?. 

In their original paper, Ginzburg and Landau noted that m certain 

materials superconductivity would not be suppressed at the thermodynamic 

critical field, H , but would persist to a field H = J 2 . KB . The 
c C2 c 

parameter n is approximately X /%, where X is the depth of penetration 

L L 

of a magnetic field m the London theory Abrikosov considered this 

phenomenon m more detail and discovered that for fields between 

H 5 H ^ //2kh =H f quantized vortices of supercurrents occurred 
ci c C C2 

m these materials The condition for a material that can have these 
vortices is that H be ^ 1/^/2 These materials are referred to as type- 1 1 
superconductors 
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Abrikosov found that m these vortices the current density vanishes 

at the center, reaches a peak value at a distance of approximately 5 

from the core, and decreased at larger radii As a consequence of this 

behavior, the magnetic field is greatest at the center of the fluxoid. 

/ - 7 2 

The flux cp m these fluxoids is h /2e or 2,07 X 10 G - cm . 

T o of 

Although the Ginzburg- Landau theory describes these properties, it 
is mathematically difficult to use The model to be presented here 
describes all these features m a mathematically tractable but approxi- 
mate form. 


Formulation of the Model 
5 

De Gennes and his coworkers developed a model of type- II super- 
conductivity starting with the London equation 


CURL(X 2 CURL ja ) + & = 0 . (l) 

For a single fluxoid with cylindrical symmetry and X = X , the London 

L 

penetration depth, equation (l) becomes 


2 d 2 b dh 

x — 2 + x T" 

dx dx 


2 , 
x h 


= 0 


( 2 ) 


m terms of x = r/X The solutions to equation (2) are the modified 

L 

Bessel functions K (x) and I q (x) which have the unphysical feature of 
a singular core A cutoff can be introduced by letting X be spatially 
dependent. This penetration depth is related to the superconducting or- 
der parameter, \jr, by X = X . /+• 

L 


Now equation (l) becomes 


2 d 2 h 


dx 


1 


3 + |l + xty 


df 


— 2 


dx 


dh 

dx 


x 2 t 2 h = 


(*') 
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If 


♦ = n-; ( 3 ) 

§ 

where y - (x + e 2 ) and N and e are field dependent variational param- 
eters (e for cutoff and N for normalization) then equation (2') becomes 


2 d 2 h dh 2 . 

y — " + y t; - y h = 
dy 2 d y 



, " \ 

Equation (2 ) has the solution 

h = C NK (Ny) + C NI (Ny) ( 4 ) 

1 0 20 

Where and C g are constants (independent of position but not field) 

By introducing a cutoff into the London equation, we have obtained 
trial functions for the order parameter, equation (3), and local magnetic 
induction, equation (4) These can be introduced into the Ginzburg- 
Landau free energy to determine N and e by a variational procedure. 

We investigated the possibility of other forms for equations (3) 
and (4) m detail, including such families of solutions as h proportional 
to an arbitrary function multiplied by a Bessel function of arbitrary 
argument, series solutions, exponentials of arbitrary functions, and so 
forth No other choice that we investigated gave tractable forms for 1)1 
and h, simultaneously 

A solution can now be constructed that approximates an Abrikosov 

fluxoid by causing the current density to vanish on a circular boundary 

of radius R This construction consists of approximating the hexagon 

on which the current actually vanishes by a circle This approximation 

6 , 7,8 

has been used previously by several authors We demand further that 

the flux contained inside the circle be the flux quantum cp Q . These 
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conditions require that dh/dx | = 0 and Jh • ds = Cp^ 

magnetic induction is given by B = cp o /rrR 2 

Using these conditions, we find 


h 

z 


- h(x) 


NCP ° Y(R,y) 

2n?L 2 e 

L 


J e 


s j(x) 


N9 c 
o 


8rr 2 X s e 

L 


♦ 


ft(R.y) 
P(R, e) 


A 


0 


= A(x) 


^ y P(R,y) + ^ 
2n\ ex P(R, e) 2nr 

Jj 


where 


7 (R,y) = I 1 (Ny R )K Q (Ny) + K a (Ny R )l 0 (Ny) 

and 

p(R,y) = I^Ny^K^Ny) - K^Ny^I^Ny) 


and 


y 


R 


(R 3 /X 2 + e 3 ) 1 


Finally, it is noted that 


J 


4ttA 8 

Li 



A 

2rrr 


Consequently, the 

(5) 

( 6 ) 

(7) 

( 8 ) 

(9) 

( 10 ) 

( 11 ) 


7 



London. Free Energy 


The usual form of the Lagrange energy density in the London theory 
is 


L = “ Ch 3 + (dh/dx) 2 ] (12) 

8TT 

Writing equation (l) in its usual form, the equation for the local mag- 
netic field is given by 


d 2 h ldh 
dx 2 xdx 


0 


(13) 


and has the solution 


h = D K (x) + D I (x) (14) 

1 o v ' 3 J v * 

At x = 0, both h and dh/dx become infinite leading to an infinite energy 
density at the core of a fluxoid The present model can be obtained from 
the formal transformation of x into Ny In this case, the Lagrange energy 
is 

L = -^ [h 2 + (dh/d(Ny)) 2 ] (15) 

which has the same form as equation (l2) At x = 0, however, the arguments 
of the K functions become Ne and the energy remains finite* The free 
energy per unit length of fluxoid is obtained by integrating the free 
energy density over the cross section of the fluxoid, resulting m the 
simple expression 



Gmzburg-Landau Free Energy 


Since a functional form of ty(x) has been assumed and a correspond- 
ing expression for h(x) has been calculated, it is possible to use these 
to obtain an approximate solution of the Gmzburg-Landau equations The 
Gmzburg-Landau free energy will be minimized with respect to the param- 
eters N and e to yield the values of these parameters which give the 
lowest energy 

The Gmzburg-Landau free energy density has the functional form 


f sh ~ f no 



tr 1 
+ + 

8tt 2m 


-ibViJi - — A i|r (17) 


from which the following two equations can be derived 


V 2 * +^2 1 


2A 3 H 2 


2 

r H 
*-■? 



Cylindrical symmetry and single- valuedness have been assumed, imply- 
ing that i|f can be written m the form - j^le 1 ^ Substituting this ^ 
with |i|f| = (Nx/y) into equation (l9), we find that equation (l9) is iden- 
tical to equation (ll) if the expressions for j and A from equations (6) 
and (7) are used 

Substituting h and ijr into equation (l7) and letting H = 

He = cp /2^/5 ttX§, and H = cp /2 tt£ 2 , the free energy m this approximation 
o o 

can be shown to be 



( 20 ) 
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It is seen immediately that the energy of the modified London model 
[equation (l5)] is only one part of the Gmzburg-Landau free energy 
The first term is the condensation energy and the last term is a part 
of the Vi[rV\jf* term. 

Integrating equation ( 20 ) over the cross section of a fluxoid and 
dividing by the cross-sectional area, the free energy per unit volume 
becomes 


8TT , v 

H 2 (F SH “ ^NO 


■t - $ 


-2N 2 (l - rr) + ~ ~ B N 3 (l - N s )An(y /«) 

il K 

C3 


+ B^l + £[ BVt i( 1+eV4 )] (21) 

C R 


The external magnetic field H = 4 tt(3/^b) (F^ - F ) can he obtained 

7 SH NO H 


from equation (2l) 


H 


N s e s 


H 


(1 - N 2 ) 


C2 


M y R A) 


H 1 

C3 


N H, 


C3 


+ “7 4 4 =‘£ 
2h y R B 




N 2 e s + 

4 /J 


H 


C3 


Ny(r, e) 

Bm^s'VbPAr, e) 2n 2 e3(R, e) 

R 


( 22 ) 


Finally, the Gibbs energy is given by 


G - G 
SH NH 


= F 


SH NO 


BH H 
4TT 8 TT 


(23) 


As B -* 0, N approaches 1 and y , 3 and y approach infinity Under 

R 

these conditions, the first and third terms of equation (22) vanish. 

The limit as B ”* 0 of the remaining two terms define the lower critical 
field H 

cl 
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(24) 


H 

Cl 


H 

C2 


~ (e 2 + l/n 2 ) 


K (0) 

0 v ' 

2k 2 sk ( e) 


2 

In thLS limit, the minimization procedure indicates that 0 approaches 
2 / k 2 yielding 


H 

cl 


H 

C2 


+ 



K q (V2/h) 
2a/5’mK i {Jz/h) 


(25) 


Results of Model 

It was found that when N and e were varied, the lowest energy at 
low fields was obtained when i}f slightly exceeded unity For this reason, 
the physically realistic condition that if be less than or equal to one 
was imposed. Some of the parameters were evaluated for several values 
of n, but detailed results are given here only for h = 5. These results 
are shown in Figures 1 through 5. 

The average Gibbs energy per unit volume of the specimen is shown 
in. Figure 1 as a function of the magnetic induction b/h . The magne - 

C2 

tization -4uM/H^ 3 was obtained numerically from equation (22) by choosing 
a value of B/H and calculating h/h The results of this calculation 
are shown m Figure 2 where they are compared with the Abrikosov magne- 
tization^ given by 


1 - h/h 

-4TTM/H = — — . (26) 

1.16(2h 3 - l) 

The difference between the local magnetic induction, h(x)Al , and 

C2 

the average induction b/h is shown in Figure 3. The local current 

C2 

density j(x) inside the fluxoid, normalized by cH /4tt\, the London 

e 
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FIGURE 2 


MAGNETIZATION, -47rm/H c2 VERSUS H7H c2 FOR a = 5 The broken line is 
the Abrikosov prediction 
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RATIO OF CURRENT DENSITY TO LONDON CRITICAL CURRENT DENSITY, 
cH c /4t rX AS A FUNCTION OF POSITION INSIDE A FLUXOID FOR VARIOUS 
VALUES OF THE MAGNETIC INDUCTION, K = 5 
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FIGURE 5 VARIATION OF THE ORDER PARAMETER, $( X ) WITH POSITION AND 
MAGNETIC INDUCTION FOR k = 5 

critical current density, is shown in Figure 4, the order parameter 
H x ) is shown m Figure 5. A list of parameters applicable to the case 
k = 5 is given m Table 1. 

Pinning of Fluxoids 

It is generally agreed that the pinning of fluxoids in supercon- 
ductors is caused by imperfections in the superconducting sample The 
variety of imperfections m metals is quite diverse, ranging from simple 
defects such as vacancy- interstitial pairs to gram boundaries, precipi- 
tates, and composition variations m alloys* Undoubtedly each of these 
contributes to some extent to the total pinning force exerted on the 
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Table 1 


PARAMETERS FOR K - 5 


B/k 
' c2 

H/fc 
' c2 

-4TTMAI 
' c2 

8rrG/k 2 

c 

2 

€ 

N 

0 005 

0 04530 

0 04030 

-0 8973 

0.082 

1.0025 

01 

04627 

03627 

-.8936 

089 

1 0055 

02 

.05104 

03104 

- . 8777 

091 

1 0113 

.05 

07389 

.02389 

- 8167 

095 

1 0292 

.1 

1192 

01924 

“ . 7206 

094 

1 0571 

2 

2152 

01525 

- 5575 

083 

1 0668 

.3 

3128 

01282 

- .4209 

.0736 

1 0378 

5 

.5088 

008849 

- 2092 

0575 

0.9126 

6 

6070 

006983 

- 1314 

.0515 

0 8239 

7 

.7051 

005151 

- .07183 

0465 

0.7157 

8 

8033 

003333 

- .03016 

0425 

0.5813 

9 

9015 

.001525 

- 006315 

0392 

0 3962 

95 

9506 

00621 

- .001049 

0377 

0 2537 


fluxoid lattice. Since so many different types of imperfections are 
present m commonly available materials, it is difficult to determine 
which imperfections contribute most to the total pinning force by mea- 
suring macroscopic quantities such as the critical current or the 
magnetization of the sample 

Experimentally, attempts have been made to reduce the number of 

types of pinning centers by carefully controlling the metallurgical 

preparation of the samples Numerous experiments of this type have 

introduced known amounts of precipitates m a host matrix, known amounts 

of radiation induced defects m well annealed samples, and known amounts 

9-20 

and sizes of grain and subgrain boundaries 
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Theoretical treatments of pinning are complicated by the inability 
to describe fluxoids that are nonsymmetric Thus, if a fluxoid is posi- 
tioned over a region of the sample m which the metallurgical, and there- 
fore the superconductive, properties vary significantly, the properties 
of the fluxoid cannot be described adequately. For this reason, valid 
calculations can only be made if the change m superconductive properties 
occurs over a distance that is large compared to the fluxoid dimensions, 
or if the change m properties occurs over a region that is small enough 
compared with the fluxoid that the fluxoid geometry can be assumed to be 
unperturbed Imperfections intermediate m size between these limits 
cannot be treated at present* These two limiting cases are considered 
in the following discussion 

Consider first the effect of a large imperfection such as a cell 
wall, gram boundary, or precipitate Such imperfections range in size 
from about 2 X 10" 6 cm to 2 X 10" 5 cm or larger depending on the purity 
and thermomechanical history of the sample (We will use the example of 
a cell wall or grain boundary to which defects and impurities have mi- 
grated ) These imperfections lead to a reduction m the electron mean- 

free-path, an increase m the resistivity, and consequently to an m- 

21 

crease m h according to the relation 

JL 

H =s H + H, = H + 7.5 X 10 3 y S p (27) 

o £ o e n 

where K is the value of H m the defect free portion of the sample, y^ 

is the electronic specific heat, and p is the normal state resistivity 

n 

Assume that through the cell wall of width 2w the value of H is 
triangular in shape, increasing from 0 to w and decreasing from w to 2w. 
Thus, dn/dx is constant from 0 to w and the pinning force associated 
with this gradient xn K is 
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f 

p _ _Sg _Bh ^ 3g Ah 
V 3h 3x ~ 3 k, W 


( 28 ) 


Letting the total length of imperfections perpendicular to the force 
be s, and the length of the fluxoid parallel to the imperfection be H, 
the volume over which the force acts is sw^ and the pinning force on the 
imperfection is 

f = 2 -!^Ahs£ (29) 

p 3 h 

where the force on each side of the wall has been included It will be 
shown later that Bg/Bh is negative for an increase m k, and therefore 
the force is attractive on both sides of the wall, drawing fluxoids into 
the wall 

Since the width of the wall does not occur m the final expression, it 
seems at first sight that only the change m energy between 0 and w or w 
and 2w is important* However, this formulation is only valid if w > 2R and 
dn/dx is small, otherwise, the energy of the fluxoid cannot be calculated. 
The radius of fluxoids is smallest at high fields and the calculation is 
most applicable At low fields, the radius increases and fluxoids are 
more likely to encounter regions of widely varying h Since a gradient 
m h is assumed to exist only m the cell wall, a force is excited and 
a current flows only in this part of the sample 

The current density m the layer w is found by equating the pinning 

force f and the Lorentz force jB/c* 

P 7 

J = Cf /b (30) 

P 


The average current density m the sample is 

j = cf sw/BA 
P 


( 31 ) 
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where A is the cross-sectional area of the sample and sw/A is the frac- 
tion of the sample m which the current actually flows ^ 


The field dependence of f arises from which was obtained 

by numerical differentiation and is shown m Figure 6 for n = 5 The 
same figure shows the result of differentiating the Abrikosov energy 
under the same conditions* The curves were not extended below 
b / h C 2 = 01 due to numerical problems with the present calculation and 
the inapplicability of the Abrikosov energy at low values of b/h 

C 2 

Neither curve goes to zero at B = 0 as might be inferred from the curves 
shown The Abrikosov result can be shown to be 



TA-8043-6 

FIGURE 6 PARTIAL DERIVATIVE OF THE AVERAGE GIBBS ENERGY WITH RESPECT TO 
K AS A FUNCTION OF THE MAGNETIC INDUCTION 
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-4 H 


(32) 


8TT SG 

H2 Sh 
c 


(1 - h/h ) 

1 C2 

(2 y? - l) 2 3 



H 

H 

C2 
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where p = 1 16. 


For large values of n, this varies approximately as 


0 - 2H/H (l - H/H ) 

8tt dG ^ _ C3_ ca 

H s Stt ” 1 16 h 


Using this approximation, the average current density m the sample 


cH 3 S 3 wMh 


3 . 32ttAH K 

c 3 



The form of this current density is shown m Figure 7 

To determine the effect of small defects on the pinning force, it 
is necessary to know the distribution of the Gibbs energy inside a 
fluxoid This is obtained using equations (20) and (23) with B replaced 
by the local induction h(x) m the latter equation* In this case, 


G - G 
SH NH 


H 


H 


^ r,d/ \ ,4/ \ /rt *i (H - h (x)) c /Ne* 53 

— [* (x) - t MA] + ~ |— 


4tt 


Hy 


(35) 


This energy density per unit volume is shown m Figure 8 as a function 
of xA/R with b/h as a parameter* It is surprising and previously un- 
recognized that although the average energy density m the fluxoid is 
negative (see Figure l), the local energy density can be positive and 
quite large* From the shape of these curves, it is seen that if a 
small normal region exists inside the fluxoid, the largest amount of 
energy will be gained when this normal region lies at the center of the 
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FIGURE 7 FORM OF THE CRITICAL CURRENT DENSITY WHEN PINNING IS CAUSED 
BY LONG RANGE VARIATIONS IN k (See equation 34) 
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FIGURE 8 VARIATION OF THE GIBBS ENERGY AS A FUNCTION OF POSITION 
INSIDE A FLUXOID FOR VARIOUS VALUES OF THE MAGNETIC 
INDUCTION, k = 5 (See equation 35) 
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fluxoid The slope of these curves gives the pinning force per unit 

volume of normal inclusion m the fluxoid as a function of position and 

magnetic induction. The maximum slope, or the maximum force per unit 

volume, does not occur at the center of the fluxoid but at some distance 

from the center To calculate the actual force, X and H must be known 

c 

We choose, for illustrative purposes, the self- cons is tent parameters 

H = 2 X 10 s Oe and \ = 7 5 X 10" 6 cm from which it follows that 
c 

§ = 1.5 X 10" s cm, H =14.6 kOe, and H = 659 Oe In this case, the 

cl 

maximum pinning force per unit volume of normal inclusion is shown m 
Figure 9 

To obtain the current density, it is necessary to know not only the 
force on a defect but the fractional volume of the sample occupied by 
these defects. For this purpose we let 

V 

f = _S (dynes/cm 3 ) (36) 

p dr V ' 


where V /v is the fraction of the volume occupied by normal inclusions, 
ir 

The current density is then j - cf^/B (statamperes/cm s ) This current 
density is shown m Figure 10 m amperes/cm 3 without correction for the 
term V /v. If m a particular sample 0.1% of the sample were occupied 
by _ normal inclusions, the current density m Figure 10 would be multi- 


plied by 10 This current density is compared with the functional 

JL 

form (l - H/H 3 )/(h/h 2 ) 8 which has been compared with experiment and 
discussed elsewhere 


Voids of the type just discussed are rare m comparison with the 
numerous other types of defects found m metals, with the exception of 
sintered materials which can contain a large number of voids It is 
of interest to determine the effect on the current density of small in- 
clusions that have a h value different than that of the host lattice 
In this case, the volume pinning force is 
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-dG/dr (1010 dynes/cm 3 ) 
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FIGURE 9 MAXIMUM PINNING FORCE PER UNIT VOLUME OF SMALL NORMAL 
INCLUSIONS INSIDE A FLUXOID FOR THE CASE k - 5, 

H c = 2 x 10 3 Oe AND X = 75 x 10" 6 cm 
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FIGURE TO CURRENT DENSITY CAUSED BY SMALL DISPERSED 

NORMAL INCLUSIONS IN A SUPERCONDUCTOR, k = 5 
Multiply values shown by volume fraction of sample contained 
in normal inclusions 




(37) 


F 

P 



V 

H 


dndr V 


A K 


where V /V is the relative volume of the material having a different 
value of h 

This derivative was evaluated numerically hut met with difficulties 
at low fields that have not been resolved at present The maximum value 
of the pinning force in the fluxoid is shown m Figure 11 The shape of 
this curve differs radically from those shown earlier A pronounced peak 
occurs at high fields It was found that although &G/dr decreased at 
high fields, the difference between values calculated for slightly dif- 
fering H values (a 0*5% change was used m obtaining the derivatives) 



FIGURE 11 VOLUME FORCE CAUSED BY SMALL DISPERSED INCLUSIONS HAVING 
VALUES OF k DIFFERING FROM THE HOST MATERIAL Multiply values 
shown by the volume fraction of the sample having the different value of k 
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increased with field It is assumed that since dG/br vanishes at H 

C3 

for the bulk of the sample, the pinning force will vanish at this field* 

It is not possible to calculate values at higher fields without 
resorting to double precision computations which time did not permit 
These calculations were considerably more difficult to handle than those 
previously discussed since the absolute energy must be known with great 
precision. The Bessel functions were calculated to an accuracy of only 
one part m 10 10 At low fields, the arguments of the functions become 
large and oscillations occurred about the asymptotic value, limiting the 
accuracy to as little as one part in 10 s At this time we are unable to 
state the cause of scatter m the results at low fields. It may be caused 
by the loss of precision, or by insufficient minimization of the energy 
with respect to N and e. 

The curve drawn m Figure 11 is a reasonable approximation to the 
calculated points and has been used to determine the average critical 
current density shown in Figure 12. This curve obviously depicts the 
well-known "peak effect” often observed m type-II superconductors. To 
obtain the actual critical current density in a sample, the current den- 
sity shown must be multiplied by the volume fraction of the sample 
occupied by the inclusions having different values of n and by Ah 

In the following section, we will show that a whole spectrum of 

values of H can occur m a sample, even m samples that are well annealed 

Consequently, xn most samples there will not be one value of Ah but 

several, and the current density will have contributions from more than 

one curve such as that shown Further, rather than having a sharply 

defined change m k there is a distribution m H about a central Ah. 

This will cause a smearing of the curve shown for most samples. The 

peak effect has been introduced into at least one niobium- zirconium 

14 

sample by radiation damage. This type of experiment introduces rather 
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FIGURE 12 CURRENT DENSITY CAUSED BY SMALL DISPERSED 
INCLUSIONS HAVING VALUES OF K DIFFERING 
FROM THE HOST MATERIAL Multiply values shown 
by the volume fraction of the sample having the different 
value of k and by A k 
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well-defined clusters of displaced atoms that should yield a fairly 
sharply defined Ah It should be noted, however, that other samples of 
different materials did not exhibit this effect 

The volume fraction of the sample occupied by these defects must 
be small for this calculation to be valid For this reason, the critical 
current density contributed by this phenomenon will be relatively small 
and the effect will be noticed only at high fields where the critical 
current density caused by the other mechanisms discussed is small 

The force exerted on the pinning site by the fluxoid is not always 
of the same sign, but can be attractive or replus ive at different positions 
For a positive An, the largest force is always attractive near the center 
of the fluxoid, but can be repulsive near the outer radius at low fields 
(this is thought to be true even though the absolute magnitude is m 
doubt for the reasons discussed earlier). At high fields, the force is 
attractive all the way to the outer radius. The form of the force at 
b/h = 0.3 is shown m Figure 13 Obviously, if there is more than 

one pinning site inside the fluxoid there is a high probability that the 
forces exerted on the two sites will tend to oppose each other and an 
average force is needed This behavior persisted in the present case 
to a field of about 0.6 H . At higher fields, the force was attractive 

C2 

all the way to the outer radius of the fluxoid. 

Conclusion 

The model developed is sufficiently accurate to permit detailed 
calculations of features of fluxoids in type-II superconductors which 
could not be calculated previously. The model can be used to determine 
all the usual physical properties of fluxoids such as energy, energy dis- 
tribution, magnetic induction, induction distribution, current density 
distribution, magnetic field, and magnetization 
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FIGURE 13 FORCE EXERTED ON A FLUXOID BY A SMALL INCLUSION HAVING A 
VALUE OF k DIFFERENT THAN THAT OF THE HOST MATERIAL FOR 
THE CASE k = 5 AND B/H c2 =03 (ARBITRARY SCALE) At high fields 
the attractive force extends to the outer radius of the fluxoid 
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This model was used to calculate the forces exerted on a fluxoid by 
several simple types of pinning centers (l) large defects with gradients 
m n, ( 2 ) extremely small voids, and (3) small inclusions with values of 
H different from that of the host lattice. These forces were related to 
the current densities that can be sustained by materials containing these 
defects* The resulting jc versus H curves are similar m form and mag- 
nitude to the curves generally observed experimentally The third type 
of pinning center indicated above had unusual properties, yielding a 

j versus H curve with a peak m j near H Further, it was found 

c C<d 

that although the strongest pinning force was always attractive for a 
positive bn in the inclusions, inclusions near the outer radius of the 
fluxoid can cause a repulsive force. 
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SPECIFIC HEAT CAPACITY 


Theoretical Analysis 

The specific heat of a homogeneous annealed sample can be charac- 
terized by the material parameters h and T =T (H=0) 

o c 

When the material is strained, we observe that the specific heat 

discontinuity m zero field is smeared out This means that the 

strained material must be characterized by a distribution m T Q 

Further, the strained material m nonzero fields displays an even 

larger temperature smearing of the specific heat than m zero field 

This means that m addition to the T distribution there is also a 

o 

distribution m h We have developed a formalism that describes and 
permits solving for such h and T Q distributions It serves also as a 
posteriori justification of the assertions 

The essential point m this theory is that the temperature varia- 
tion of the specific heat near the transition is extremely sensitive 
to variations in material parameters Away from the transition, the 
specific heat is particularly insensitive and, indeed, not predictable 
from theory with enough accuracy to determine material parameter dis- 
tributions 


To obtain a useful specific heat expression which displays the dis- 
continuity, the normal electronic specific heat, which at low temperature 

equals y T, is subtracted from the total specific heat C* This function 
e 

for particular values of T , H = H (T = o), H is 

o' o c 


C - 


f(T, H, T , H , h) = 


y T 
e 


f (T, 
1 


H, T 


H 


h) Ml 


(T^h) - T) 


(38) 
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where C is the specific heat The Heaviside unit function M(t (h) - t) 

c 

expresses the discontinuity The functional form away from the discon- 
tinuity is contained m f One of the better expressions for f is 

25 

that of Van Vijfeijken 


* x (t, H > V V h) = 


K 


H - Log H 


o/ T \ H 

3(— ) + 

T o a/5kH 


- 1 


(39) 


The measured specific heat m a strained material at fixed H would 

include contributions from all values of T and n present m the material, 

o 


f(T, H) = J* dTdHP T (T Q )P H (H)f i (T,H, T q , H q , h) M(T o (H)-T). (40) 


The experimental data imply that distributions m two parameters 

are required If the microscopic basis for the T and K distributions 

o 

were understood, it might result that these distributions were related 

The sensitivity of the discontinuity to the material parameters 
is best identified m the derivative df/dT, which would be infinite for 
a homogeneous, annealed material 


df 

dT 



~ df l 

f dT dn P (t ) P (h) 

J O T o K S 

0 

it (T - t) - f 6 (T - t) 

dT v c 7 1 c 


(41) 


In the region where df/dT is large, only the second term makes a non- 
negligible contribution, 


— jw - fdT dn P (T ) P (h) f, 6 (t - T) 

dT J o T o' h v ' l v c 


( 42 ) 
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Further reduction of equation (42) requires a functional form for 
Tc(H, T , h) For most materials, m the vicinity of T q 


H 





(43) 


At the transition H = H equation (43) inverts 

ca 


T (H, T , k ) = T 
CO o 



(44) 


This form applies to a variety of type-II superconductors Use of more 
precise and more complicated forms for equations, (43) and (44) has only 
a minor effect on the results below 

Further reduction of equation (42) must be performed before data can 
be analyzed In zero field, the argument of the 6 function has only a 
T q dependence The temperature dependence becomes 

« P (t) f dHP (n)f (T, k) 

T o * 1 

H = 0 

In nearly every form conceivable for f x , the variable T Q would 

occur as T/T . Therefore equation (45) is P (T) multiplied by a con- 
o T 

o 

stant that is later absorbed m the normal izat i on # Since we want P (t) 

0 

normalized such that its mtergal is unity, 

J' dT o P T (T o ) =1 (46) 

o 

equation (45) is an algorithm for determining P (t) directly from the 

T o 

experimental data 


(45) 


df 

dT 


35 



When the magnetic field is not zero, a modification of the above 

approach must be used to determine the H distribution Note that m 

equation (44) for T , K appears only m the form Ha/ 2H 0 = . This 

implies that the observable is really H , hence our analysis of the H 

C2 

distributions will be m the form of H distributions Equation (42) 


is changed to 


c 2 


df 

dT 


- /dT dH P m (T ) P_ (H ) f 6(T - T) 

J o Cg T q v o H c 2 1 c 

C 2 


(42 0 


and equation (44) becomes 


T C (H, T 0 , 


3_ 



(44') 


First examination of equation (42 0 would suggest that we have a hopeless 

unfolding problem, because elimination of the one 6 function will still 

result in an integral over P and P Thus, we would have at best, 

T H 

0 C 2 

a Fredholm integral equation of the first kind for P 

H 

C 3 

This is not the case, however, for elimination of the delta func- 
tion m equation (42 0 gives us an integral over two distributions, each 
of which is fairly narrow or peaked To examine their relative contri- 
butions to the derivative of the specific heat given by equation (42 > ), 
we can take the derivative with respect to temperature of the argument 
of the function as given by equation (44 0 The result is 


AT 

T 



H/H AH 

X Cg ££ 

2 1 - H/H H 

c 2 eg 


(47) 
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The important point is that the contribution to the width of the tran- 
sition consists of two parts, one due to the distribution m T and one 

o 

due to the distribution m H The latter distribution has its width 

c 2 

multiplied by a factor that is quite large at high fields For example, 

at field values near the largest value for which measurements were taken, 

2,700 gauss, the contribution of the T o distribution to the width of the 

specific heat transition was less than 5$ of the contribution due to 

the H distribution. This means that m the evaluation of equation 
c 2 

(42 ) (for the purpose of analysis of the data m high fields), it is 

quite satisfactory to treat the T q distribution as a delta function, 

P = 6(T -T ) We then obtain the following equation for the H dis- 
T o ° ° c 2 

tnbution 


df 


dT 


2 



(48) 


where 


H * = 



(49) 


Equation (48) and (49), 


plus the normalization condition P dH = 1, 

J H c 2 

c 3 


provide a complete algorithm for analyzing specific heat data for H 

C2 

distributions. Because of the considerations implied by equation (47), 
it is worth while to do such an analysis only at large values of fields 
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Specific Heat Measurements, Data, and Analysis 


The heat capacity of high purity niobium samples was measured using 

the apparatus shown m Figure 14 This device uses an ac calorimetry 

26 

technique developed by Sullivan and Seidel 

In the apparatus shown m Figure 14, a laser beam which is chopped 
at a known frequency, U), is used to heat a sample The temperature of 
the sample rises above that of the copper reference block by an amount, 

AT , determined by the equation at the bottom of the diagram. The 
equation states that this temperature rise consists of two parts ( l) 
a dc temperature rise given by Q^/k and, (2) an ac temperature rise pro- 
portional to l/c. This temperature change, At, is measured by using a 
gold- iron versus copper thermocouple. The thermocouple is mounted on the 
sample and coupled to a copper reference block whose temperature is mea- 
sured using a dc thermocouple connected to a Keithley- type 147 nano volt- 
meter. The actual increase in dc temperature above the temperature of 
the block is determined from the thermocouple by using the nanovol tmeter 
An ac signal is produced m the thermocouples by the temperature rise 
and is amplified through a low temperature transformer The signal is 
measured by a PAR HR-8 lock- in amplifier tuned to the same frequency, w, 
as the chopped laser beam 

It should be noted that there are some correction terms m the ac 
temperature rise that do not make it strictly proportional to l/C. By 
judicious choice of 0), t (relaxation time to block determined by the 
thermal conductivity of the thermocouple wires), and (thermal relaxa- 
tion time across the sample itself), the correction terms can usually 

be made less than 1$. Thus by properly choosing cd , T , and T , the ac 

1 2 

temperature rise can be measured accurately However, to make an abso- 
lute measurement of the heat capacity, Q q must be known accurately 
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An effective measure of the Q is made by using a standard sample 

o 

of known heat capacity (usually copper) m one of the two chambers shown 

in Figure 14 The laser light, and its associated Q , produce an ac 

o 

signal This signal is assumed to be inversely proportional to the heat 

capacity of the known sample By measuring the ratio of the two ac sig- 

« 

nals with the same Q^, and at the same temperature, an absolute measure 
of the heat capacity of the unknown sample can be determined 
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THERMAL POST BLOCK WITH HEATER 
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TA-8043-38 

FIGURE 14 SPECIFIC HEAT APPARATUS 

As primary advantages, this technique has the ability to measure 
relative changes m heat capacity as a function of temperature with a 
precision of about 0 1$, and the ability to make continuous measurements 
of heat capacity as the temperature is varied thereby increasing the 
speed of the data-takmg process Furthermore, since the mass of the 
additional apparatus required to sense the temperature and heat the 
sample has been minimized, an absolute measurement of the heat capacity 
of samples whose weights are approximately 150 mg or less, can be made 
with an accuracy of approximately 3$ 
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The temperature resolution is approximately 1 mK while the absolute 


accuracy is approximately 20 mK The accuracy of the magnetic field 
measurement is approximately 1 % which implies a 28 gauss probable error 
at the highest field measured of 2,745 gauss 

The niobium specific heat capacity data is shown m Figures 15, 16, 
17, and 18 as measured for four different field ranges near zero, 500, 
1,300, and 2,700 gauss Four different samples were measured. The first 
was an annealed sample. The second had a strain of 0.07. The third had 
a strain of 0.20, and the fourth had a strain of 0.75. The strains were 
introduced by rolling The measurements m zero field show a systematic 
increase of transition temperature with strain, and a systematic increase 
of width of the specific heat transition with strain The results near 
500 gauss resemble the results of zero field with the exception of 
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FIGURE 15 SPECIFIC HEAT OF STRAINED AND UNSTRAINED NIOBIUM SAMPLES 
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the shifting of the relative positions of the curves due to the varia- 
tions of magnetic field strength The change of character that is noted 
in the results at 1,300 and 2,700 gauss is due to the shift of the tran- 

i 

sitxon width from a transition that is dominated by the T distribution 

o 

to one dominated by the H or H distribution as described m equation (47) 

C2 

T Distribution 
o 

Taking the derivative of the specific heat with respect to the 

temperature from the data shown in Figure 15 for zero fields yields the 

T distribution as given m equation (45) 
o 

The T distribution of the annealed sample is a delta function or 
o 

a gaussian distribution with a standard deviation of less than 0 5 mK 

Figure 19 shows the T 0 distribution for a strain of 0 07. This distnbu 

tion is fit reasonably well by assuming a gaussian distribution with 

T =9.225 K and a standard deviation of 4 mK When the strain reaches 0.20, 
o 

the distribution is broadened as shown in Figure 20. This distribution is 

still sufficiently well fit by a gaussian distribution with T 0 = 9 268 K 

and a standard deviation of 11 mK Note that there is a high temperature 

tail that extends several standard deviations above T . When the strain 

o 

is as large as 0 75 as shown m Figure 21, we found that the distribution 

of T was broadened considerably and not even close to a gaussian dis- 
o 

tribution m shape The distribution extends many times as far above 

the mode of the transition temperature, T =9 285, as it does below T 

o 0 

The distribution to higher temperatures appears more exponent ial~l ike 

than gaussian The behavior of the modes, T q , versus strain as shown in 

Figure 22, indicates that for small strain there is a rapid increase m 

T with strain, and for large strain there is an asymptotic approach to 
o 

a maximum value at e = 1 # 0 It is useful also to show how the width of 
the transition changes with strain m Figure 23 The width is defined 
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FIGURE 19 DISTRIBUTION IN T Q FOR A NIOBIUM SAMPLE 
OBTAINED FROM SPECIFIC HEAT CAPACITY 
MEASUREMENTS IN ZERO FIELD The strain 
was 0 07 
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FIGURE 21 DISTRIBUTION IN T Q OBTAINED FROM SPECIFIC HEAT CAPACITY 
MEASUREMENTS IN 2ERO FIELD The strain was 0 75 
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FIGURE 22 CHANGE OF MODE TRANSITION TEMPERATURE 
WITH STRAIN IN NIOBIUM 
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FIGURE 23 CHANGE OF WIDTH OF THE T Q 
DISTRIBUTION WITH STRAIN IN 
NIOBIUM 
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as the half width at a relative height of 0 135 of maximum height For 
a gaussian distribution such a definition of width would yield twice the 
standard deviation The width increases rapidly with strain for small 
strain, but shows a much lower derivative for large strain 

H or k Distributions 
C3 

Analysis of the high field specific heat measurements near 2,700 
gauss, according to the algorithm of equations (48) and (49), yield the 
H or K distributions shown m Figures 24, 25, 26, and 27 Figure 24 

C2 

shows the n distribution for zero strain The results were surprising 

because they indicate two narrow, clearly separated distributions m H 

The relative separation is approximately 2$ From earlier work with a 

27 

less quantitative model, we had expected such a behavior to be more 
likely m highly strained samples m which there is a well-defined cell 

wall structure The lower peak with H = 3,333 has approximately 33% 

C3 

of the area of the entire H _ distribution. The upper peak occurs at 

c 3 

H =3,395 The K scale uses H -1,993 gauss for niobium as obtained 

C« 0 

28 

by Finnemore, Stromberg, and Swenson 

In the sample with 0 07 strain shown in Figure 25, the same basic 
structure persists except that each of the peaks have been broadened 
and the relative magnitude has changed For 0 07 strain H = 3,293 

Cjg 

and 3,358 We do not regard the shifts of peak position as significant 
because a 37-gauss shift, which is only slightly more than the probable 
error for the field measurements, would cause the two structures to be 
m agreement 

The sample which was strained to 0.20 and is shown in Figure 26, 
demonstrates that the distributions have been broadened further into 
two plateaus and the structure at the lower value of H was nearly de- 
pleted at the expense of the peak at the higher value of K Figure 27 

for a strain of 0 75 shows a behavior m the H or h distribution that 

C2 
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FIGURE 24 DISTRIBUTION IN H c2 OR k IN AN 

ANNEALED NIOBIUM SAMPLE OBTAINED 
FROM SPECIFIC HEAT CAPACITY 
MEASUREMENTS IN HIGH FIELD The 
field was 2670 gauss 
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FIGURE 25 DISTRIBUTION IN H c2 OR k IN A NIOBIUM 
SAMPLE OBTAINED FROM SPECIFIC HEAT 
CAPACITY MEASUREMENTS IN HIGH FIELD 
The field was 2745 gauss The strain was 0 07 
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FIGURE 26 DISTRIBUTION IN H c2 OR k OBTAINED FROM SPECIFIC 
HEAT CAPACITY MEASUREMENTS AT HIGH FIELD The 
field was 2745 gauss The strain was 0 20 
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NIOBIUM 
H « 2745 
e = 0 




FIGURE 27 DISTRIBUTION IN k OBTAINED FROM SPECIFIC HEAT 
CAPACITY MEASUREMENTS AT HIGH FIELD The 
field was 2745 gauss The strain was 0 75 




is less complicated in structure The peak that remains corresponds to 

the highest H or h peak m the annealed sample Further, the lower K 
C3 

material has disappeared and has been converted into a high H distribution 

that extends about 5$ above the value of H The extension of the dis- 
cs 

tribution at high strain toward a higher K clearly represents improvement 
of the magnetic field properties of the material, as well as improvement 
m the pinning properties as discussed earlier 


Conclusion 

Our major conclusion from the specific heat capacity data analysis 
is that the technique developed for the analysis of high resolution spe- 
cific heat data has implications to the analysis and understanding of the 
metallurgy of superconducting materials that extend far beyond the par- 
ticular samples analyzed in this report There appears to be an oppor- 
tunity for a detailed understanding of the metallurgy, and a characteri- 
zation of the contribution of different types of imperfections and defects 
to probability distributions over the relevant parameters We have ob- 
served m detail the tendency of imperfections to increase the transition 
temperature and H or H We have observed also the strain induced 

C2 

broadening of the distributions m T and H, causing a cutoff below 

0 

the unstrained values of T and H and extending the tail of the dis- 

o C2 

tributions some considerable distance above T and H These distribu- 

o c 2 

tions clearly indicated how large strains can improve the temperature and 
magnetic field characteristics of superconductive materials 

This analysis technique should be applied to practical materials 
such as superconducting magnet wire at various stages of manufacture to 
determine the effects of the manufacturing processes on the distributions 
of and n Because our time and funds were depleted sooner than ex- 
pected, we were unable to perform extensive metallurgical studies on the 
niobium samples This would have enabled us to identify the types of 
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defects in the sample, and determine methods of improving the T and 

o 

H distributions However, for practical assessment of materials we 

C2 

have provided a basis for completely characterizing superconducting 
materials by macroscopic specific heat measurements 

In fact, the technique appears useful enough that T and H or 

o C3 

H distributions should be treated as primary data for understanding the 
behavior of superconducting materials 

Finally, the and H distributions can be interpreted as the 
relative fraction of volume that has a particular T q or H What we do 
not know is whether these changes m properties are concentrated m 
small regions of the sample or whether they are finely distributed 
throughout the sample m small clusters. This question must be answered 
by detailed microscopic metallurgical examination 
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BITTER PATTERN MEASUREMENTS 


The purpose of the Bitter pattern experiments was to observe the 
distortion of the fluxoid lattice m imperfect type-II superconducting 
niobium resulting from the interaction of the fluxoxd lattice with the 
defect structure of the material. For example, if the fluxoid lattice 
maintained its hexagonal symmetry, but showed a continuous variation m 
lattice parameter or fluxoid spacing, this could be related to a flux 
gradient Such flux gradients result from pinning, and direct observation 
on the scale attainable with this technique would allow correlation of 
these flux gradients with microscopic defect distributions The surface 
structure of the sample observed is of primary importance m that varia- 
tions m the microscopic sample thickness would introduce flux gradients 
resulting from variation of the total free energy of the system with the 
sample thickness These observations were to act as a guide to calcula- 
tions made using the theoretical model developed m this work 

Experimental Method 

The method of observing the distribution of flux m type-I and 

type-II superconductors is a modified Bitter pattern technique The 

29 

Bitter pattern technique was developed to allow the observation of 
magnetic domain walls m iron, nickel, and cobalt In this technique, 

o 

a colloidal suspension of small ( less than 200-A diameter) ferromagnetic 
particles is placed on the surface of a ferromagnetic material The 
particles experience a force when m the presence of an inhomogeneous 
magnetic field and are drawn to the points on the surface of the sample 
at which the magnetic field intensity is highest In this colloidal 
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suspension, the particles are free to move and distribute themselves m 
a manner determined by the field intensity distribution at the sample 
surface By observing the variations in the density of the particles, 
the spatial variation of the magnetic field (i e , the position of the 
domain walls) can be determined. 


This method has been applied also to investigations of the inter- 
mediate state m type~I superconductors, using both ferromagnetic and 

superconducting powders This application has been reviewed by Livingston 
30 

and De Sorbo In these observations, the normal or the superconducting 
state domains are delineated by the magnetic powders More recently a 

31 

technique with which higher resolutions are possible has been reported 
This technique is similar to one developed independently at SRI , It con- 
sists of evaporating a ferromagnetic metal (iron, cobalt, or nickel) from 
a refractory heater m a helium gas atmosphere onto a sample m the super- 
conducting state The ambient pressure is extremely important and de- 
termines the resulting particle size Fine particles with diameters down 

o 

to 80 A can be produced m this manner at gas pressures between 0 5 and 

1 0 torr These particles are allowed to settle on the superconducting 

surface and distribute themselves m a pattern determined by the local 

magnetic flux density. Individual fluxoids as well as the fluxoid lattice 

32 33 

m type-II materials have been observed by this technique. ’ 


The resolution necessary for the observation of fluxoids is defined 
by the value of the parameter K from the GLAG theory and, therefore, by 
the penetration depth and the coherence length The spacing, d, of these 
fluxoids m the hexagonal array is described in terms of the bulk magnetic 
induction, B, and is given by l/d 2 = qB/cp 0 , where q is a geometric factor 
describing the packing density m the hexagonal array, and cp Q is the flux 
contained m a single fluxoid For relatively pure niobium, having a 
H of 1.2 at 4 2 K, a fluxoid spacing of the order of 1,700 A is expected 
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for inductions of approximately 700 G This represents about 20 magnetic 

particle diameters and allows sufficient resolution for observation of 

individual fluxoids This spacing increases as the magnetic induction, 

B, decreases and larger spacmgs are possible with materials of higher 

H due to the resulting lower values of H In fact, spacmgs of the 

Cl 

o 

order of 6,300 A are predicted for inductions of the order of 50 G. 

In the Bitter technique, particle agglomeration is necessary for 

the observation of the distribution of the magnetic flux m the sample 

34 

The conditions for this to occur have been determined by Kittel for 
the case of a colloidal suspension of particles He found that this dis- 
tribution was field-dependent and given as P(h) = sinh ( (iH/kT)/(p,H/kT) , 
where P(h) is the particle density expected for a field H, \i is the mag- 
netic moment of a particle, k is Boltzmann T s constant, and T is the ab- 
solute temperature Analysis of this equation shows that usable particle 
densities are expected when ( pH/kT) £ 3 0 For iron particles of 80-A 
diameter at a temperature of 4°K, this requires the local field to be 
larger than about 4 G This condition is met since the field at the 
fluxoid core is greater than the applied field. 

Experimental Apparatus and Procedures 

The experimental apparatus was designed so that samples at tempera- 
tures between 2.0 and 5 0°K could be studied. The apparatus generally 
consisted of a vacuum-tight thxn-wall stainless tubular probe 24 m. 
long X l-l/4 m diam that was isolated from the liquid helium bath 
The sample was mounted in this pi'obe The probe was filled to one at- 
mosphere pressure with helium gas after pumping and purging three times. 
It was then inserted into the liquid helium bath after cooling by immer- 
sion in liquid nitrogen. The probe fit (with 0 020-m clearance) into 
a superconducting niobium wire solenoid (458 gauss/A) 
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A schematic drawing of the inserted probe is shown m Figure 28 
The sample (typically a disc 0 250-m diam and 0 005 to 0 020-m. thick) 
is mounted at ( l) and is thermally grounded to an OFHC copper sample 
mounting block (7)* This block can be removed from the end of the probe 
for mounting and examination of the sample A vacuum tight seal is made 
using a sheared indium gasket (8) The sample is optically shielded (2) 
from the evaporation source (3). The evaporation source (3) is mounted 
on the end of two 0.150-diam stainless steel rods (4) that act as power 
leads for the evaporation source and are isolated from the probe walls 
(6) by teflon spacers ( 12) The sample is positioned at the center of 
the magnet. The whole system is immersed m liquid helium contained m 
a double-wall glass dewar ( 10) . 

A typical experimental sequence is as follows The sample is posi- 
tioned as described above The helium gas pressure m the probe is re- 
duced to 250 (i This initial pressure was determined empirically and is 

o 

the pressure necessary for particles of approximately 200 A to be formed 

during the evaporation. The sample is allowed to sit in a constant mag- 

netic field for one hour to ensure complete deposition of the ferromagnetic 

powder so that no erroneous patterns will be observed because of powder 

deposition during removal of the probe from the liquid helium. The sample 

is then removed from the liquid helium and allowed to warm to room temper- 
ed 

ature over a period of three hours The sample is removed and a 500-A 
carbon layer is deposited. Subsequent observations are made on samples 
m this configuration 

Experimental Results 

A total of 19 experiments were performed, eleven of which are listed 
m Table 2 The remaining eight are not described m detail since they 
were directed toward empirical definition of the experimental conditions 
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(1} SAMPLE 

(2) OPTICAL HEAT SHIELD 

(3) COBALT WRAPPED 
TUNGSTEN FILAMENT 
EVAPORATION SOURCE 

{41 STAINLESS STEEL POWER 
LEADS AND EVAPORATOR 
SUPPORT 

(5) NIOBIUM WIRE SUPER- 
CONDUCTING MAGNET 
(458 gauss/amp) 

(6) STAINLESSSTEELTHIN 
WALL EVAPORATION 
CHAMBER 

(7| ANNEALED 0 F H C 
THERMAL GROUND 
SAMPLE MOUNTING 
BLOCK 

(81 INDIUM GASKET 

(9) LIQUID HELIUM 

(10) DOUBLE WALL 
GLASS DEWAR 

(11) MAGNET SUPPORT 
STRUCTURE 

(12) TEFLON CENTERING 
PIECES 


TA-8043-10 


FIGURE 28 


EXPERIMENTAL APPARATUS FOR BITTER 


PATTERN MEASUREMENTS 
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Table 2 


BITTER PATTERN EXPERIMENTS PERFORMED ON DEFORMED AND ANNEALED NIOBIUM 


a> 

to 


Deposition Conditions 


Experiment 

Temperature 

Field 

Time 

Current 

Pressure 


Observations 

Number 

<°K) 

( gauss) 

( sec) 

(A) 

(n) 

Optical 

Electron Microscope 

8043-1 

to 

o 

1,000 

0 4 

55 

200- 

X 

X 

e=0 





1,500 

rv 


8043-2 

4 2 

1,000- 

4 

55 

200- 

X 

X 

o 

CM 

O 

II 

<£> 


500 



1,100 



8043-3 

4 2 

1,000- 

4 

60 

200- 

X 

x - 

e = 0 20 


500 



350 



8043-4 

4 2 

1,500- 

4 

60 

215- 

X 

X 

e = 0 75 


500 



2,000 



8043-5 

4 2 

1,000- 

4 

60 

220- 

X 

X 

e = 0 07 


500 



700 



8043-6 

4 2 

1,000- 

4 

60 

230- 

X 

X 

0 = 0 07 


800 



750 



8043-7 

4 2 

1,000- 

5 

60 

225- 

X 

X 

6 = 0 20 


800 



400 



8043-8 

4 2 

1,000- 

5 

60 

225- 

X 

X 

o 

El 

Ci> 


500 



750 



8043-9 

4 2 

1,000- 

5 

60 

200- 

X 

X 

CO 

II 

O 

-0 

U1 


500 



400 



8043-10 

4 2 

750 

7 

60 

200- 

X 

X 

0 = 0 





700 



8043-11 

4 2 

1,000- 

7 

60 

180- 

X 

X 

0 = 0 75 


0 



350 





necessary to produce the correct helium gas pressure m the sample probe 
and of the evaporation conditions necessary to produce the correct amount 
of powder without unnecessary heating of the samples. 

Samples deformed from 0 to 75% by rolling were examined Typically, 
the external field was raised to 800 to 1,000 gauss and then reduced to 
approximately 500 to 800 gauss. This sequence was used so that stable 
flux distributions were established m the samples and so that the trapped 
flux density m the samples exceeded the external field. The powder is 
attracted to the sample and controlled experiments are possible. 

All samples listed m Table 2 were examined using both electron 
and optical microscope techniques, e.g , during optical microscop y, obser- 
vations were made of sample 8043-5, a pure niobium sample deformed 7% by 
rolling. A steady applied field of 500 gauss was maintained on the 
sample after cycling the field to 1,000 gauss Since the lower critical 
field of high purity niobium at 4 2°K is nearly 1,500 gauss, the sample 
was never fully in the type- IX state and the mixed state structure was 
defined poorly. The structure of the sample is shown m the low magni- 
fication micrograph of Figure 29 Essentially, no microscopic structure 
is obtainable m this picture 

In Figure 30, the flux distribution at the sample edge is shown. 

A quasi-mtermediate state structure is apparent m this figure with the 
regions of high flux density lying nearly parallel to the sample radius. 
This structure is consistent with flux distributions obtained m other 
type-II materials at fields less than the lower critical field 

In Figures 31, 32, and 33, the flux distributions near the center 
of the sample are shown. Two specific characteristics are apparent. 

First, the flux distribution is nonuniform, but locally isotropic What 
appear to be flux spots of nearly uniform size are shown. These consist 
of high densities of individual quantum fluxoids. The local field 
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FIGURE 29 MICROGRAPH OF SAMPLE STRUCTURE 
(SAMPLE 8043-5) 
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FIGURE 30 FLUX DISTRIBUTION AT THE SAMPLE 
EDGE. Lines of high flux density show 
a quasi-intermediate state structure. 
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FIGURE 31 


FLUX SPOTS AND REGIONS OF RATHER HIGH 
FLUX DENSITY NEAR THE CENTER OF THE 
SAMPLE 
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FIGURE 32 


FLUX DISTRIBUTION AT THE SAMPLE 
CENTER 



FIGURE 33 FLUX DISTRIBUTION AT CENTER OF SAMPLE. 

The regions of high flux density have a structure 
that is not fully resolvable in this micrograph. 
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intensity m these flux spots is of the order of 700 gauss, the equili- 
brium induction for fields just above the lower critical field of mo- 

o 

bium Therefore, the fluxoid spacing expected is of the order of 2,000 A 
which is beyond the resolution of optical microscopes. This is shown m 
Figure 33 where obvious, though unresol vable , structure is seen within 

o 

the individual flux spots At a magnification of 500X, 2,000 A is equal 
to 0.1 mm Therefore, it is likely that individual fluxoids may be shown 
m this micrograph because of some randomness m the powder deposition, 
but the detailed structure cannot be resolved. 

The carbon film was stripped from this sample using standard repli- 
cation techniques and observed m the electron microscope in transmission 
Essentially none of the deposited powder was removed and, therefore, no 
data were obtained This problem was encountered in all the samples 
investigated, and although the scant optical observations indicated that 
interesting patterns were formed, we were unable to observe them m 
detail. 

The experimental problem of removing the deposited powder from the 
sample surface for observation m the electron microscope limited the 
observations m this work. Several techniques typically used to aid m 
the removal of carbon films from sample surfaces were tried These in- 
cluded etching, electropolishmg, and deposition of a parting layer before 
the formation of the powder pattern. None of these was successful, al- 
though the inclusion of a parting layer did show indications of leading 
to a solution to the problem 

Conclusion 

A series of experiments directed toward the direct observation of 
flux distributions m type-II superconductors by a modified Bitter pat- 
tern technique have been conducted The purpose of these experiments 
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was to observe quantum fluxoid distributions m imperfect type-II super- 
conductors and to correlate imperfections m the fluxoid lattice with 
defect structure 

Although optical micrographs indicated that several interesting 
Bitter patterns were formed on the niobium samples, the inability to 
strip the powder from the sample surface prevented us from observing 
these patterns with the electron microscope We were unable to resolve 
this problem with the time and funds available 
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CONCLUSIONS 


• The theoretical model formulated m this project is a good 
approximation to the Gmzburg-Landau theory and is mathe- 
matically tractable 

• Flux pinning arises from a variety of metallic defects, each 
of which contributes to the critical current density m a 
different way. Three types of pinning were analyzed 

• The peak effect was found to be caused by small dispersed 
inclusions of material with a k value higher than that of 
the host material 

• A method for analysis of the h and T 0 distributions m super- 
conductors caused by deformation was formulated The method 
uses experimental specific heat measurements* Using this 
analysis, distributions of h and T 0 m annealed and deformed 
niobium samples were determined 

• Bitter patterns showing the location of fluxoids m pure 
annealed and pure deformed niobium were formed. Problems 
associated with replicating these patterns precluded 
observing the patterns m detail using electron microscopic 
techniques 
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RECOMMENDATIONS 


• The theoretical model developed in this program provides a 
simple, accurate means of calculating the properties of 
type-II superconductors and should be explored m more detail 
The results of the flux pinning calculations should be com- 
pared with experiments on carefully controlled samples The 
peak effect analysis should be checked experimentally by 
introducing a microscopic dispersion of a second material 
with a k slightly higher than that of the host material 

If the predictions can be verified, a method of controlling 
the current density at high fields would result 

• The technique for analysis of specific heat data developed 
in this program should 'be applied to superconducting wire 
and other practical materials The technique can aid 
considerably m the understanding of the relationship between 
the macroscopic performance and the microscopic metallurgy 

of superconducting materials 

• The Bitter pattern technique used shows promise of providing 
a direct measure of the pinning forces exerted on fluxoids 

by different types of pinning centers. This information would 
be valuable in assessing the accuracy of theoretical models 
of the type developed m this program This technique should 
be developed to its full potential. 
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